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ABSTRACT 

We construct a simple and robust approach for deriving constraints on magnetic fields in galaxy 
clusters from rotation measure (RM) maps. Relaxing the commonly used assumptions of a correlation 
between the magnetic field strength and the plasma density and of a power-law (in wave number) 
magnetic field power spectrum, and using an efiicient numerical analysis method, wc test the consis- 
tency of a wide range of magnetic field models with RM maps of 11 extended sources in 5 clusters, 
for which the data were made available to us. Wc show that the data reveal no indication for a radial 
dependence of the average magnetic field strength, and in particular no indication for a correlation 
between the gas density and the field strength. The RM maps of a considerable fraction of the sources 
either require or are consistent with the presence of a spatially uniform magnetic field of a relatively 
small strength, 0.02-0.3 /iC, which contributes significantly to the RM. The RM maps of all but one 
source do not require a power-law magnetic field power spectrum, and most are consistent with a 
power spectrum dominated by a single wave length. The uncertainties in the magnetic field strengths 
(and spatial correlation lengths) derived from RM maps exceed an order of magnitude (and often 
more). These uncertainties imply, in particular, that there is no indication in current RM data for 
a systematic difference between the magnetic field strengths in radio-halo clusters and in radio-quiet 
clusters. With the improvement expected in the near future of the quality and quantity of RM data, 
our analysis method will enable one to derive more accurate constraints on magnetic fields in galaxy 
clusters. 
Subject headings: cxtragalactic magnetic fields, galaxy clusters 



I. INTRODUCTION 

The strength and spatial structure of magnetic fields 
in the intra-cluster medium (ICM) of galaxy clusters 
have been estimated using a wide range of methods. 
These methods include analyses of the properties of ra- 
dio relic and halo sources, of cluster non-thermal X-ray 
emission, of cluster cold fronts (appearing in X-ray im- 
ages), and of Faraday rotation measure maps of polar- 
ized radio sources (for reviews, see Carilli & Taylor 2002; 
Govoni & Feretti 2004). The results of these analyses im- 
ply that the ICM of all galaxy clusters are permeated by 
~ 1 /xG magnetic fields. Such fields play a critical role 
in determining the energy balance in the ICM plasma 
through their effects on heat conduction and cosmic ray 
propagation, with implications to many cluster phenom- 
ena (for reviews, see Feretti et al. 2012; Ferrari et al. 
2008). 

In this paper we discuss the constraints that are im- 
plied by Faraday rotation measure (RM) data on the 
properties of the ICM magnetic fields. The Faraday ef- 
fect causes a rotation of the plane of polarization of an 
electromagnetic wave, which is linearly proportional to 
the component of the magnetic field in the direction of 
propagation, chosen henceforth as z. The RM is defined 
as 
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where ip is the angle of rotation, A is the (radio) wave- 
length, RMq = 0.812 radrn"^. rig is the (free) electron 
number density, and the integral is along the line of sight 
to the source (Burn 1966). Cluster RM maps are ob- 
tained by fitting the wavelength dependence of the po- 
larization angles measured along each line of sight to the 
A^ lawof Eq. (1). 

The two dimensional RM maps do not allow one to 
uniquely determine the three dimensional magnetic field 
structure. Moreover, RM maps are available over only 
a small fraction of the angular extent of the clusters. 
Thus, RM constraints on the ICM magnetic fields are 
typically derived under simplifying assumptions regrad- 
ing the magnetic field structure. Commonly used as- 
sumptions are that the magnetic field strength is corre- 
lated with the plasma density, and that the power spec- 
trum of the spatial distribution of the magnetic field is 
a power-law in wave numbers. Furthermore, the range 
of model parameters (e.g. the range of wave numbers 
or the range of power-law indices) explored is limited in 
most derivations of magnetic field constraints, either by 
assumption or by some qualitative analysis of the two 
dimensional RM maps, in order to reduce the required 
computational resources. 

The analysis presented here of the constraints imposed 
on the ICM magnetic fields improves on earlier work in 
several ways. First, we relax the commonly used as- 
sumptions of a correlation between the magnetic field 
strength and the plasma density, and of a power-law (in 
wave number) power spectrum. Indeed, we find that the 
RM maps show no evidence for a correlation between 
the field strength and the plasma density (as previously 
shown, e.g., by Eilek & Owen 2002), and that a power- 



law power spectrum is generally not require to account 
for the observations. Second, we use a more efficient nu- 
merical analysis method, that allows us to explore a wide 
range of model parameters. Finally, we apply a uniform 
analysis to sources observed in several clusters, which al- 
lows us to draw conclusions applicable to the population 
of clusters (rather than to a single cluster or source). 

RM maps have been derived and analyzed for radio 
sources in 16 nearby galaxy clusters based on observa- 
tions of more than 40 radio sources (see Table 1). The 
RM data available and earlier analyses of this data are 
described in § 2. Our analysis method is described in § 3, 
and our results are described in § 4. We have applied our 
analysis method to all the data that were made available 
to us, 11 sources within 5 clusters, by F. Govoni, A. 
Bonafede, V. Vacca, and their collaborators. We thank 
them for making the data available, and note that our 
analysis method may be straightforwardly applied to all 
existing data. Our results are summarized and their im- 
plications are discussed in § 5. In particular, we address 
the issue of magnetic field strength bi-modality between 
merging and non-merging galaxy clusters. 

In our analysis we assume that the observed RM are 
entirely due to the ICM magnetic fields and not intrinsic 
to the sources (e.g. Pizzo et al. 2011), and that the con- 
tribution of the magnetic field of our own Galaxy to the 
RM is smaller than the fit errors (e.g. Taylor & Perley 
1993). We adopt a ACDM cosmology' with Hq = 
72kms-iMpc"\ ^u = 0.27, and r^A = 0.73. 

2. EXISTING DATA & EARLIER ANALYSES 
2.1. Existing RM maps 

RM maps have been derived and analyzed for radio 
sources in 16 nearby galaxy clusters (Table 1). De- 
tailed RM maps were constructed for more than 40 radio 
sources (Table 2) at different impact parameters (angular 
distance between the line of sight and the cluster X-ray 
center), lying within or behind the ICM. As can be seen 
in Figure 1, the radio sources observed in these studies 
are distributed over a wide range of distances from their 
clusters' centers, with roughly half of the sources within 
2rc {re as defined by Eq. (5)). 

The sources included in our study are indicated in Ta- 
ble 2. More than half (7/11) of our analyzed sources are 
within Coma, a prominent radio halo (RH) cluster, and 
4 of these sources are coincident with the RH. We also 
have one source in each of the clusters A2065, A2142, 
A2199 and Ophiuchus, none of which are hosting a RH. 

2.2. Earlier analyses 

For five of the galaxy clusters listed in Table 1 
(Coma, A2199, Hydra A, A2255, and A2382), de- 
tailed models of the ICM magnetic fields have 
been constructed and constrained by the observed 
RM maps (Taylor & Perley 1993; Govoni et al. 2006; 
Guidetti et al. 2008; Bonafede et al. 2010; Vacca et al. 
2012). For the rest of the sources only a limited anal- 
ysis was carried out (Bonafede et al. 2011; Govoni et al. 
2010). We briefly summarize below the main results of 
both types of studies. 

2.2.1. Clusters for which detailed analyses were carried out 



For Coma, Bonafede et al. (2010) assumed that the 
spatial distribution of the magnetic field is given by a 
realization of a Gaussian random field with a Fourier 
power spectrum that is a power-law of the wave num- 
bers. They have assumed a Kolmogorov spectral index, 
n = 11/3, and by comparing various statistical properties 
of the observed RM images and those simulated, the spa- 
tial wavelength range was constrained to A,„i,i = 2 kpc 
and A,nax = 34 kpc. 

In order to generate a correlation of the field strength 
with the plasma density, the (real space) magnetic field 
obtained in the realizations was multiplied by an r- 
dependent factor to generate 



(S(r))=BoK(r)/no]^ 



(2) 



where the average is over all realizations. The RM maps 
derived from the magnetic field realizations were com- 
pared to the observed RM maps, to constrain the central 
magnetic field strength Bq to be in the range 3.9-5.4 /iG, 
and the correlation index rj to be in the range 0.4-0.7 

(iCTC.L). 

A similar analysis has been carried out for A2199 by 
Vacca et al. (2012). In their analysis, Amax = 35 kpc, 
Amin = 0.7 kpc, and n = 2.8 were used, based on fitting 
a power-law spectrum to the 2D RM map (from which 
Amax = 35 ± 28 kpc, Aniin = 0.7 ± 0.1 kpc, and n ~ 
2.8 ± 1.3 were inferred). Using these values and following 
a method similar to that of Bonafede et al. (2010) for 
Coma, they find Bq ^ 11.7 ± 9.0 /xG and r/ = 0.9 ± 0.5 
(Icrc.L). 

Besides Coma and Abeh 2199, the ICM mag- 
netic fields of only three other clusters - Hydra A 
(Taylor & Perley 1993), A2255 (Govoni et al. 2006), and 
A2382 (Guidetti et al. 2008) - have been investigated 
with a similar approach, i.e. by simulating random mag- 
netic field models and constraining their properties by 
comparison to the RM data. 

For Hydra A, several models were explored for the ICM 
fields. Taylor & Perley (1993) find - 20 ^iC to - 40 ^iC 
magnetic field strength assuming spatial variability on 
r^ 4 kpc scale, and some evidence for magnetic fields with 
cluster scale ordering. 

For AbeU 2255, a RH cluster, Govoni et al. (2006) es- 
timate the cluster core field strength to be ~ 2.5 /iG. 
They try several models with a few different values for 
the power-spectrum power-law index, keeping the field 
strength at the cluster center constant at 2 — 2.5 ^G, and 
conclude that it is necessary to use a (spatially) variable 
power index (in the range n = 2 — 4). 

For Abell 2382, Guidetti et al. (2008) produce a fit 
with central magnetic field strength, So ~ 1-13 /iG [la 
confidence region) for models with different correlation 
index between the average magnetic field strength and 
the gas density. Note that all their models produce sim- 
ilar magnetic field strengths when averaged over a large 
volume (1 Mpc' ). 

2.2.2. Clusters with limited analyses 

Comparing the standard deviation of each RM source, 
crm (averaging the observed values over all valid pixels 
in each source), to the distance from the cluster center, 
not normalized to account for the different cluster sizes 
(e.g. to Tc, the core radius of the /3-model), and assuming 
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t Mass ref.: Girardi et al. 1998, /3-model ref.: Govoni et al. 2001. 

J Mass ref.: Demarco et al. 2003, /3-model ref.: Guidetti et al. 2008. 

Tab. 1. — Galaxy Clusters used in RM studies, including some relevant cluster properties. (1) Cluster name; (2) Redshift (Nasa/ipac 
EXTRAGALACTIO DATABASE); (3) Merger activity (Bonafede et al. 2011, and references therein); (4) Diffuse radio emission: G = giant 
halo, M = mini halo, R = relic (Ferctti et al. 2012); (5) Cool-Core: W = weak, S = strong (Hudson ct al. 2010); (6) Virial mass 
(Rciprich & Bohringer 2002, unless otherwrise noted); (7) (8) (9) /3-modcl (Eq. 5) parameters with 1-a fit errors adapted according to our 
adopted cosmology (Chen et al. 2007, unless otherwise noted). 
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Fig. 1. — Distribution of RM sources. Each marker represents the mean position of an extended RM source, relative to its host galaxy 
cluster (in the sense that the radio emission passes through its ICM) X-ray center. Note that a few of the central sources cover a relatively 
large area of their respective cores. 



a "Gaussian cell modeP", Govoni ct al. (2010) find these 
galaxy clusters are consistent with Bq ^ 1-4 /iG. 

3. METHOD 

3 an analytical formulation based on the approximation that the 
magnetic field strength is constant throughout the cluster and that 
its direction varies over a single length scale 



We assume that the cluster's magnetic field may be 
described as a realization of a random process, the sta- 
tistical properties of which we seek to constrain using 
the measured RM map. We consider a statistical model 
acceptable if the power spectrum of the observed (two 
dimensional) RM map is consistent with being drawn 
from the distribution of power spectra predicted by the 
model. To determine whether or not this is the case, we 



TABLE 2 
Radio sources analyzed in previous work 
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t The radio souree is coincident with the cluster RH. 

J The radio source might be coincident with the cluster RH (edge of diffuse emission). 

Tab. 2. — The radio sources for which RM maps have been derived and analyzed in previous works, along with some of their relevant 
characteristics. In this work we analyzed the sources in Coma, A2065, A2142, A2199, and Ophiuchus (italicized). (1) Host GC name; 
(2) Source name; (3) Radio morphology: (H)T = (Head) Tail, (N/W)AT = (Narrow/Wide) Angle Tail; (4) Fanaroff-Riley classification; 
(5) Source object type; (6) Impact parameter relative to X-ray center; (7) Line of sight position: oo = background source, = cluster 
member; (8) RM reference: 1. Bonafede et al. 2010; 2. Eilek & Owen 2002; 3. Feretti et al. 1999; 4. Govoni et al. 2001; 5. Govoni et al. 
2006; 6. Govoni et al. 2010; 7. Guidetti et al. 2008; 8. Perley & Taylor 1991; 9. Pizzo et al. 2011; 10. Taylor &; Parley 1993; 11. Taylor et al. 
2001; 12. Taylor et al. 2002; 13. Vacca et al. 2010; 14. Vacca et al. 2012. 



compare the distribution of some measure of the differ- 
ences between the RM power spectra of different real- 
izations of the model with the distribution of the same 
measure of differences between the model's power spec- 
tra and the observed power spectrum. We define the RM 
power spectrum as 



RM(fc)' = (|RM(fc„fc,y)r 



I [kx .ky) 



(3) 



where {kx,ky) is the 2D Fourier fc-vector corresponding 
to the projected distance on the plane of the sky, and 
the average on the rhs is over all vectors with magnitude 
k. For our analysis we use two measures of the distance 
between two power spectra, RMi(fc) and RM2(fc): the 
sum and the maximum of the absolute value of the log- 



arithm of the ratio of powers in the two power spectra 
at all wave numbers, i.e. X^fc |lc)g(RMi(fc)/RM2(fc))| and 
maxfc |log(RMi(fc)/RM2(fc)j|, respectively. For each sta- 
tistical model we generated r^ 500 realizations and com- 
pared the resulting distance measure distributions using 
the Kolmogorov-Smirnov (KS) test. Similar results are 
obtained for the two distance measures used. 

We write the magnetic field as a sum of Fourier com- 
ponents, ^j^ B(k) exp(ik • x) with 

B(k) = Ske*^f^ ici cos(afj) 4- ka sin(af2) , (4) 

where ki and ka are two unit vectors orthogonal to k 
(chosen to align along the z-direction of Cartesian coor- 
dinates) and to each other (ensuring V • B = 0), and 
i?k is real. Since the magnetic field is assumed to be 



statistically homogeneous and isotropic, which choose 
the phases a^ and P^ to be uniformly distributed in 
[— TT, +7r]. The fc-spacc resolution is chosen as ~ l-KJrc- 

We consider two types of power spectra for the am- 
plitude JSj.: a single characteristic wave length A, with 
i?k = for fc 7^ 27r/A, and a power-law spectrum with 



index n, with k^B\ 



for fc in some fixed range 



[fcmin,fcmax] (corresponding to A,„in = lkpc,An 
32 kpc) and zero elsewhere. In both cases, we allow the 
presence of a uniform (fc = 0) component (as proposed 
by Taylor & Perlcy 1993), with amplitude i?honi- The 
average magnetic field strength, Bq = (3 (i?f ))^'^, is as- 
sumed to be constant throughout the cluster (i.e., with 
no radial dependance). As described in § 4, the simple 
models we consider are sufficient to account for the vast 
majority of observations, and hence more general models 
are not motivated by the data. 

For any realization of the magnetic field model, the 
Faraday RM map is obtained using Eq. (1) and a j3- 
model density profile. 



= 7io(l + rVr2) 



-ip 



(5) 



from the literature (as listed in Table 1). We treat clus- 
ter members as situated precisely halfway through the 
ICM, and other sources as lying at infinity. For these 
two cases, the integral of Eq. (1) may be solved analyt- 
ically, as shown in Appendix A. Note also that due to 
the linear nature of Eq. 1, it is possible to explore as 
many spectral index values as desired, without having to 
simulate the magnetic field more than once at each wave 
number k. 

In order to account for the finite angular resolution of 
the observations, RM values at each pixel are obtained 
by properly averaging the values obtained in a band of 
surrounding pixels corresponding to the reported angu- 
lar beam size. Our results are not sensitive to ~ 10% 
changes in the beam size adopted for the calculation. 
RM values are derived in the simulations only at coordi- 
nates specified by pixels with observed RM data, and at 
a band of pixels surrounding them. It is thus not neces- 
sary to use a large cubic Fourier grid and not necessary 
to use heavily resource consuming FFT calculations. 

Finally, after the convolution with the beam, Gaus- 
sian noise is added, with standard deviation equal to the 
mean observed RM fit errors. Our realizations showed 
that the noise has only a little influence on the RM 
power-spectrum. Some contribution to the power on 
short wavelengths is obtained only when the level of the 
noise is increased to ^ 100%. 

4. RESULTS 

In this section we present the results of our analysis, 
which was described in § 3. For every source we sim- 
ulated five types of magnetic field models: single-scale, 
power-law, both with and without a homogeneous com- 
ponent, and a simple homogeneous field. Figures 2 and 3 
show the RM power spectra of 5C4.114 and A2199A, and 
compare them to those obtained for two types of mod- 
els: a homogeneous field model for 5C4.114 and a single 
scale model for A2199A. The results are summarized in 
Table 3, where 95% c.l. ranges arc shown for accepted 
model parameters (as is the convention also in what fol- 
lows). 



For some of the sources, most notably 5C4.114 and 
A2142A, the RM power-spectra resemble closely those 
obtained for a homogeneous magnetic field, reflecting the 
fact that the power spectrum is determined by the shape 
of the sources. Our realizations show that a homoge- 
neous sub /iG magnetic flcld is sufficient in order to gen- 
erate their observed RM power-spectra: 0.026-0.030 /iG 
for 5C4.114 and 0.18-0.20 /uG for A2142A. 

Figure 4 shows the minimal magnitude of the mag- 
netic field implied by observations of each source, which 
is the minimal magnitude of a homogeneous field required 
to generate the maximal (absolute) RM value of each 
source. Using Eqs. (1) and (5), this minimal field is 



Bq > CLosVomax 



|RM(x,y)| 1 



/3- 



\P) 



RMonor,V^r(|/3 - i) 

(6) 
where clos ~ 1 for a background source or clos = 2 for a 
cluster meinber, and (a;, y) correspond to right ascension 
and declination relative to the cluster center. The lower 
limits are are ~ 30% higher than the acceptable ranges of 
values inferred above for 5C4.114 and A2142A. The smafl 
offset supports the validity of our simulations' results. 
It is due to the fact that the lower limits are derived 
based on the RM of a single pixel, in contrast with the 
(weighted) averaging over many pixels inherent to our 
statistical analysis. 

For most sources, a single-scale model is consistent 
with the data. All accepted ranges of field strengths of 
our single-scale models, with and without a homogeneous 
component, are plotted in Figure 5 as a function of the 
projected distance from the cluster center. We also indi- 
cate 

BcMB = (STraTcMB) ^ 3 /iG (7) 



as a dashed gray line, where relevant, and notice that 
accepted field strengths are distributed around this line. 



TABLE 3 

Simulation results - accepted models 



(1) 


(2) 


(3) 


(4) 


(5) 


Source 


Comp. 


-Bhom 


Bo 


A 






[mG] 


[mG] 


[kpc] 


5C4.042 


U-l-P 


0.005-0.025 


3.3-7.6 


_t 


5C4.074 


S/P 


< 0.57 


0.3-19 


1.3-32 


5C4.081 


S / P 


< 0.28 


0.52-5.9 


2-32 


5C4.085 


S/P 


< 0.46 


0.52-7.7 


3.3-32 


5C4.114 


U / S /P 


< 0.085 


< 1.7 


1-32 


5C4.127 


U+S / P 


< 0.41 


1.0-19 


1.5-32 


5C4.152 


S/P 


< 1.3 


0.17-16 


1-32 


A2065A 




no models were 


: accepted 




A2142A 


U / S /P 


<0.3 


0.2-3700 


3.5-32 


3C 338 


S/P 


< 0.23 


0.45-3.9 


11-28 


OPHIB 


S/P 


< 0.19 


0.5-24 


4.9-32 



t Acceptable power-law indices for the source 5C4.042 range 
from n = to n = 2.4. 

Tab. 3. — (1) Source name; (2) Acceptable models (95% c.l.): 
U=uniform, S=single-scale, P=power-law, '/'=or, '-|-'=and. Note 
that a uniform component may always be added to a varying 
one; (3) Uniform component magnetic field strength; (4) Average 
magnitude of the spatially varying magnetic field component; 
(5) Allowed spatial variation scales of the single-scale model. 
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Fig. 2. — A homogeneous field model for 5C4.114. Left panel: a histogram of realization power-spectra, for a homogeneous field model 
with -Bhom = 26.1 nG, compared to the data power-spectrum (red). Colors indicate the number of realizations for which the various values 
were obtained. Right panels: the (binned) CDF and PDF (top and bottom respectively) of the sum of absolute logarithmic differences 
between the data to the model (red) and between realizations (black). The maximized KS test statistic is Pks = 0.016. 
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Fig. 3. — Same as Figure 2, but for the single-scale field model with A = 16 kpc and Bq 
a homogeneous component). The maximized KS test statistic is pxs = 0.762. 



((3^2^)1/2 = 0.964 AtG for A2199A (without 



All the radio sources analyzed in this work, except 
A2065A (the only source for which no model was ac- 
cepted and may thus require a more complicated model), 
are consistent with some power-law magnetic field in the 
range of scales of Amin = 1 kpc to A^ax = 32 kpc, chosen 
to accommodate scales down to the angular resolution 
and up to the typical source size. In Figure 6 we show 
the ranges of acceptable magnetic field strengths, for our 
power-spectrum models, along with the analytical lower 
bounds found using Eq. (6). Again, there is no apparent 
radial trend of magnetic field strength, and the values 
are scattered around the -Bcmb fine. 

We note that although there is a correlation between 
the spectral index of the (3D) magnetic field power- 
spectrum and the resulting shape of the (2D) RM power- 
spectrum (e.g. the average RM power spectrum of a 
single-scale model peaks at the same scale as that of the 
magnetic field power spectrum) , the spectral index of the 
magnetic field power spectrum is not given by a power- 
law fit to the RM power spectrum. 

5. DISCUSSION 

We presented a simple and robust approach for de- 
riving constraints on magnetic fields in galaxy clusters 
from rotation measure (RM) maps (§ 3). Relaxing the 
commonly used assumptions of a correlation between 
the magnetic field strength and the plasma density (e.g. 
Bonafcde et al. 2010; Vacca et al. 2012) and of a power- 
law (in wave number) magnetic field power spectrum, 
and using an efficient numerical analysis method, we 
tested the consistency of a wide range of magnetic field 
models with RM maps of 11 extended sources in 5 clus- 
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Fig. 4. — The lower-bounds on the magnetic field, as derived 
using Eq. (6) (arrows), and the ranges of magnetic field strengths 
in the homogeneous field models accounting for the A2142A and 
5C4.114 RM data (95% c.l. error bars). The small offset between 
the bounds and allowed ranges is due to the fact that the lower 
bounds are derived based on the RM of a single pixel, in contrast 
with the (weighted) averaging over many pixels inherent to our 
statistical analysis. 



ters (Coma, A2065, A2142, A2199, and Ophiuchus), for 
which the data were made available to us (§ 4). This con- 
stitutes a significant extension of previous work (§ 2.2). 

We have shown that the data reveal no indication for a 
radial dependence of the average magnetic field strength 
(see fig. 5), and in particular no indication for a correla- 
tion between the gas density and the field strength. The 
RM maps of a considerable fraction of the sources either 
require or are consistent with the presence of a spatially 
uniform magnetic field of a relatively small strength (see 
table 3), 0.02-0.3 ^G, which contributes significantly to 
the RM. The presence of such spatially uniform fields 
has already been suggested by Taylor & Perley (1993) in 
their analysis of Hydra A, but was not taken into ac- 
count in recent studies (Taylor & Perley 1993, infer field 
values which are ~10 times larger than our estimate of 
the uniform components in Coma, but the RM values in 
Hydra A are also an order of magnitude larger) . 

The RM maps of all but one source do not require a 
power-law magnetic field power spectrum, and most are 
consistent with a power spectrum dominated by a single 
wave length (see table 3). The magnetic field values in- 
ferred by earlier analyses (§ 2.2) are consistent with our 
results. However, we find that the uncertainties in the 
magnetic field strengths (and spatial correlation lengths) 
derived from RM maps are significantly larger than those 
inferred by earlier analyses, which considered a more lim- 
ited range of models and model parameters. We find that 
the uncertainties exceed an order of magnitude and of- 
ten more (see table 3). This implies, in particular, that 
there is no indication in current RM data for a system- 
atic difference between the magnetic field strengths in 
radio-halo clusters and in radio-quiet clusters. We note 
that the spatial characteristics of the 2D RM maps are 
highly affected by the beam resolution and by the shape 
of the underlying extended sources. 

Our results have important implications also to the 
discussion of the mechanism responsible for the ob- 
served correlation (Brunctti et al. 2007) between the ra- 
dio power and the thermal X-ray luminosity of radio 
emitting galaxy clusters, which is a subject of consider- 
able debate. Several models for the diffuse synchrotron 
emission, the so called radio halos (RHs), have been sug- 
gested, with different assumptions regarding the origin 
of the emitting electrons. In some models the emit- 
ting electrons (and positrons) are secondaries generated 
by p-p interactions between the CR protons and ther- 
mal ICM protons (e.g. Dcnnison 1980; Kushnir et al. 
2009; Kcshct & Loeb 2010), whereas other models as- 
sume that the emitting electrons are re-accelerated by 
turbulence from a preexisting population of non-thermal 
seed electrons in the ICM (secondary or otherwise, e.g. 
Brunctti et al. 2001; Petrosian 2001). As pointed out by 
Kushnir et al. (2009), the upper envelope of the corre- 
lation between radio-power and X-ray luminosity is nat- 
urally obtained in models where the synchrotron radi- 
ation is produced by secondary electrons, provided the 
magnetic field strength exceeds i?cMB = 3/iG. In such 
models, the suppression observed in radio quiet galaxy 
clusters of the synchrotron luminosity, by a factor of ^ 10 
compared to the upper envelope of the correlation, must 
be due to lower values of the magnetic field, B < 1 ^G. 
Therefore, accurate measurements of the magnetic fields 
in radio emitting and radio quiet galaxy clusters may dis- 
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Fig. 5. — The accepted ranges of magnetic field strengths of the varying component, Bq = {3(^3'^'))^''^, for single-scale models. The 
top-left and top-right panels depict the accepted ranges of magnetic field strengths as a function of the impact parameter for the single-scale 
models without and with a homogeneous component respectively. Note that the y-axis scale differs between panels. In the bottom-left and 
bottom-right panels wc show the corresponding BqA^^^ values. The error bars represent the range of accepted values (95% c.l.), and the 
arrows indicate the lower-bounds as derived in eq. 6. 



criminate between models (see, however, Basu 2012, for 
lack of such bi-modality between RH and radio quited 
clusters in Planck data). 

The range of acceptable magnetic field strengths de- 
rived in this work implies that it is difficult to test the 
hypothesis of magnetic field bi-modality using RM maps 
(see the scatter around B = Bcmb in figs. 5 & 6). Note, 
in particular, that we have shown here that there is no 
evidence to support some of the assumptions used in ear- 



lier works (discussed in § 2.2, e.g., Bonafede et al. 2011; 
Govoni ct al. 2010), which argued against the possibility 
of such magnetic field strength bi-modality. 

A further comment is in place here regarding the possi- 
ble magnetic field bimodality. Based on a statistical anal- 
ysis of the radial trend of depolarization, Bonafede et al. 
(2011) find no evidence for a difference between popu- 
lations of galaxy clusters with and without a RH. They 
find, using the logrank test, that the null hypothesis that 
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Fig. 6. — The accepted ranges of magnetic field strengths of the varying component, Bo = (3 (Bf ))^''^, for power-law models. The left 
and right panels depict the accepted ranges of magnetic field strengths for power-law models without and with a homogeneous component, 
respectively. Note that the y-axis scale differs between the panels. The error bars represent the ranges of accepted values (95% c.l.), and 
the arrows indicate the lower-bounds as derived in Eq. (6). 



the values of depolarization for both populations (with 
and without a RH) are different realizations of the same 
distribution is acceptable. However, one should note that 
their analysis is based on RM measurements of sources 
that lie at distances from the cluster center that extend 
to '^ lOrc, while radio halos typically extend only up to 
~ 2rc- It is unclear what the conclusion of the analysis 
would be if only sources up to this shorter distance are 
considered. 

Accurate measurements of the magnetic fields in 
galaxy clusters might be possible in the near future, 
when the quality and quantity of RM data is improved. 
Higher sensitivity, resolution, and imaging capability will 
be provided by the EVLA extension of the VLA. Most 
important for RM maps, the resolution will be improved 
to 0.02-0.2 arcseconds in the cm wavelength range, al- 
lowing one to explore the contribution of shorter spatial 
fluctuation scales. LOFAR will explore the universe at 
radio-frequencies below 250 MHz, and will thus be able 
to detect diffuse cluster radio sources which are brighter 
at lower frequencies, in particular radio halos with very 
steep spectra (Feretti ct al. 2012). The longer wave- 
lengths also make it possible to detect lower RM values, 
allowing one to fill some of the missing areas in current 
RM maps. 

New wide-field surveys of the sky will offer one the op- 
portunity to carry out RM studies of many more galaxy 
clusters. Two such projects are planned with the ASKAP 
telescope. EMU is a radio sky survey project which will 
make a deep (^ 10 ^Jy r.m.s.) radio continuum survey 
covering the entire Southern Sky. It will have about 45 
times the sensitivity of the NVSS, and an angular res- 
olution 4.5 times better (Cassano et al. 2012). Because 
of the excellent short-spacing UV coverage of ASKAP, 



EMU will also have higher sensitivity to extended struc- 
tures such as cluster halos (Cassano et al. 2012). In POS- 
SUM, the plan is to use ASKAP's unique survey capa- 
bilities to measure the Faraday rotation of three million 
cxtragalactic radio sources over 30,000 square degrees 
("POSSUM Wide"). 

APERTIF, the new Phased Array Feed system that 
will be installed on WSRT, will increase by about a factor 
30 the observed area on the sky, at frequencies of 1.0 GHz 
to 1.7 GHz (Cassano ct al. 2012; Rottgering et al. 2011). 
The survey speed of APERTIF, and many of the other 
characteristics, will be very similar to ASKAP. The ex- 
tremely large field of view of APERTIF would enable the 
WODAN project, which has been proposed with the aim 
of charting the entire accessible northern sky at 1.4 GHz 
down to 10 /iJy r.m.s. and '^ 1000 dcg down to 5 /iJy 
(Cassano et al. 2012). 

Other instruments, such as the LWA, MeerKAT, and 
the African SKA, will also contribute to increasing the 
quantity and quality of RM observations. 

As RM data becomes more abundant, with more 
sources being observed, closer to and within galaxy clus- 
ter cores, our proposed analysis method will provide 
higher quality estimates of the magnetic field strength. 
Furthermore, it may prove useful to exploit the extra in- 
formation encoded in the polarization maps themselves 
(from which the RM maps are derived), as has recently 
been done for A2199 (Vacca ct al. 2012), in order to bet- 
ter constrain magnetic field models. 
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APPENDIX 
A. FROM ICM MAGNETIC FIELDS TO RM MAPS 
We derive below an analytic expression for the RM, eq. (1), obtained for a magnetie field of the form given by eq. (4), 



B(r) = NAkAuA(t)Bo ^ k^'^Fk 



with 



+ 1 +7r 



Fk= J2 J2 '^"'" kiCos(«fJ+k2sin(aj-j) 



^ikr 



(Al) 



(A2) 



Here, A0, Aw, Afc are the fc-space resolution of the azimuthal, (cosine of the) polar, and radial (spherical) coordinates, 
respectively, and A^ is a normalization factor to be determined below by requiring Bq = (3 (i?f ))^/^. Writing 






+ik-r 



u—0 {p— — 7r 



jp- 



-ki cos(a_k) — k2 sin(Q;_k) ) e 



-ik-r 



and choosing a_^2 = '^n ^i^d P_f-^ + /^^j = tt, so that B(r) is real, we have 

+ 1 +7r 



EE(i 



ki cos(aA) + k2 sm(a/ 



U—0 — — TT 
+ 1 +7r 



=+<^n+k-r) ^ g-</3^+k-r) 



2 E E (^1 cos(af2) + ka sin(ajj) cos(/3j-^ + k • r) 



It— (j) — — -K 



and 



B^(r) = 2NBo ^ cos(af2) cos(/3fj + k • r) /s^'^AfcAuA^Vl - u^. 



(A3) 



(A4) 



(A5) 



where the factor vT— ?? arises from the z component of ki (the z-componcnt of ka vanishes). The average magnetic 
field strength is then given by 

{B,{v)f) = AN^Bl J2 Ma^] cos(a,-,J) (cos(/3,-,^ +ki • r) cos(/3^^ +k2 • r 



ki,k2^ 



- 2 "ki,k2 



=hs. 



2<'ki,k2 



Jl fc, ^Afc,AM,A0,y'l-w| 

'i=l,2 

= ^2^2 ^ A:^-" (AfcAMA(/))' (1 - u2) , 

which is independent on r. Defining A^ = N' / (AkAuAcj)), 

{Bl) = N'Bl Y, k^-''AkAuA(j) (l - u^) 

pk„,^„ pl 

I d2 / lA-n 



N'B, 



femin 



/fc'*-"dA: / (l-u^)du 







27r 



1,5-n _ jL5-n o 

Af'S2:^Hl2^ 5HiiL£27r 



we find 



AT' 



5 — n 3 
1 5-n 

«'max fi-tnin 



(A6) 



(A7) 
(A8) 
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and the final expression for the magnetic field is 



B4r)^Bo^ J2 J2 J2 NkNnVl - u^ cos(aj^) cos(/3j^ + k • r) 



k=k„,i„ u=0 4>=0 



Here 



Ni 



(5 - n) fc4-"Afc 

L.5-n _ 1,5-n 
"-max n'.^jr. 



N^ ^ AuAcI). 



The rotation measure, RM, is now given by 

RM{x,y) = RMo / B.Uedz 



k,u,(p 



v2 — / dz eos(/3j^ + X + kuz) I 1 



x'^ + y^ + z"^ 



-6 



=/ 



where 

One may write 

where 



b~—(3 , X = 'sv 1 — M^ (a;cos(/) + y sin( 



lo = (r/rc) 



/ = /o3? 



1-3/J 



V 277 J 



C = z/?* , w = kur , r = yr^ + a;^ + ; 



For a background source, the integral over the line of sight ( J_ ) is a Fourier transform. 



-F, 



{(i + C^)-''}H = ^M^-^/^,_.(M), 



(A9) 

(AlO) 
(All) 

(A12) 

(A13) 
(AM) 
(A15) 

(A16) 



where K is the modified Bessel function of the second kind, and the final expression for the rotation measure is 

'^kNn \. . — 



RM(a;, y) = RMoBoUor, ^ N^Nn ^^3^^ ( ^ ) i^3^_ 1 (fcur) 



X cos 



r(i/3) V 2r ; 

(ttfj) cos( /3j=j + fcv 1 — u^ (xcos0 + ysin0) j 



(A17) 



For a cluster member the integral over the line of sight (/„ ) is a Fourier Cosine/Sine transform. 
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l-b 



,6-:! 



2r(5) 



K^_l{uj)+iTT 



2sin(67r) 



(A18) 



where / is the modified Bessel function of the first kind, and L is the modified Struve function, and the final expression 
for the rotation measure of a cluster member is 



RM(x, y) = RM„B„nor, ^ N^N, 



Vl-u^ fk 



k.u,<f) 



r(|/3) V 2r 



X <^ cos(/?jj + x) Kip_i {kur) + sin(/3j^ + x) tt 



2sin(|/37r) 



(A19) 
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